Hierarchical Fractal Weyl Laws for Chaotic Resonance States in Open Mixed Systems 
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In open chaotic systems the number of long-lived resonance states obeys a fractal Weyl law, which depends 
on the fractal dimension of the chaotic saddle. We study the generic case of a mixed phase space with regular 
and chaotic dynamics. We find a hierarchy of fractal Weyl laws, one for each region of the hierarchical decom- 
position of the chaotic phase-space component. This is based on our observation of hierarchical resonance states 
localizing on these regions. Numerically this is verified for the standard map and a hierarchical model system. 
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It is just a century ago that Hermann Weyl published his cel- 
ebrated theorem on the asymptotic distribution of eigenmodes 
of the Helmholtz equation in a bounded domain fl] which has 
found fundamental applications in the context of acoustics, 
optical cavities and quantum billiards ||2|-|4J]. For a quantum 
billiard with a d dimensional phase space the number A/'(/c) of 
eigenmodes with a wave number below k is on average and in 
the limit of large k given by JV{k) ^ /c^/^ up to corrections of 
higher order p-9]. Only recently, this has been addressed for 
open scattering systems, where for the case of fully chaotic 
systems a fractal Weyl law was found |IiqI42Q|] . Due to the 
opening of the system one classically obtains a fractal chaotic 
saddle (sometimes also called repeller), which is the invariant 
set of points in phase space that do not escape, neither in the 
future nor in the past [21, 22T. Its fractal dimension S plays 
an important role quantum mechanically: The number A/" of 
long-lived resonance states is given by a fractal Weyl law. 
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which here is stated for open chaotic maps, where the k de- 
pendence is replaced by the dependence on the effective size 
of Planck's cell h. 

Generic Hamiltonian systems exhibit a mixed phase space 
where regular and chaotic motion coexist [23 J. Regular reso- 
nance states of the open system obey a standard Weyl law, and 
naively one would expect that the chaotic resonance states fol- 
low the fractal Weyl law, Eq. ([T]). This ignores, however, that 
the dynamics in the chaotic region of generic two-dimensional 
maps is dominated by partial transport barriers, see Fig. [T] 
They are hierarchically organized with decreasing fluxes to- 
wards the regular regions [|24N28l] . They strongly impact the 
system's classical 1124143 (Jl and quantum mechanical [3ll44in 
properties, and lead to, e.g., the localization of eigenstates in 
phase space Oil 
tions ll34l35[ l39l 



33, 37, 4in and fractal conductance fluctua- 



Classically, the chaotic saddle in generic systems gives rise 
to an individual fractal dimension for each region of the hier- 
archical decomposition of phase space (29] . These are effec- 
tive fractal dimensions as they are scale dependent and coin- 
cide in the limit of arbitrarily small scales. Quantum mechan- 
ically, fractal Weyl laws for open systems with a mixed phase 



space have been investigated in Refs. ll42l445ll . but the influ- 
ence of the hierarchical phase- space structure remains open. 
In particular, the individual effective fractal dimensions of the 
chaotic saddle have not been taken into account, so far. 

In this paper we propose a generalization of the Weyl law to 
open systems with a mixed phase space. We find hierarchical 
fractal Weyl laws, 
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one for each phase- space region Aj of the hierarchical de- 
composition of the chaotic component in a generic two- 
dimensional phase space. Here, 5j denotes the effective frac- 
tal dimension of the chaotic saddle in each region. Quantum 



_._,Ailv— ^M|»^^ 




10'' 



10^ h ^ 10^ 

^-2 



FIG. 1. (Color online) Hierarchical fractal Weyl laws Kfj vs. h~^ 
(filled symbols) counting hierarchical resonance states in the outer 
(Ao, triangles) and inner (Ai, circles) chaotic regions, enclosed by 
partial barriers (solid lines, inset), for the standard map dX i^ — 2.9, 
and corresponding typical Husimi representations (right) for h — 
1/1000. They are compared to the box-counting scaling N^ vs. £~^ 
(open symbols) giving distinct fractal dimensions 5j of the corre- 
sponding regions of the chaotic saddle (inset). The rescaled count- 
ing functions are given by Nj{h) — Nj(h) ■ h ■ fj and N^^{s) — 
Nf^e) • e^, respectively. 



mechanically, this result is based on our observation of hierar- 
chical resonance states, which predominantly localize on one 
of the regions Aj . Their number J\fj follows the hierarchical 
fractal Weyl laws, Eq. (O. This is confirmed for the generic 
standard map and a hierarchical model system. 

Classical properties: We first review the classical proper- 
ties of the chaotic saddle in a generic mixed system and illus- 
trate them for the prototypical example of the Chirikov stan- 
dard map [46] . It is obtained from the kicked rotor Hamilto- 
nidin H{q^p^t) = T{p)-\-V{q) XlnGZ^(^~^) with kinetic en- 
ergy T{p) = p'^/2 and kick potential V{q) = -^ cos(27rg'). 
At integer times t it leads to the symmetrized map g^+i = 
^t +TXp*),Pt+i = P* - IV'iQt^i) withp* =pt- lV\qt) 
on the torus [0,1) x [— ^,^). We open the system by defin- 
ing absorbing stripes of width 0.05 on the left and right, see 
Fig. \T\ This leads to a chaotic saddle F, for which a finite 
time approximation is shown in Fig. [T] for Hi = 2.9. The regu- 
lar islands (including small embedded chaotic layers), which 
appear on all scales in phase space, form an invariant set of 
positive measure. The chaotic saddle F of the open system is 
strongly structured by the presence of partial barriers. They 
originate from Cantori or stable/unstable manifolds of hyper- 
bolic periodic orbits [28]. Partial barriers provide a hierarchi- 
cal tree-like decomposition ll27ll of the chaotic component of 
phase space into regions Aj. Figure [T] shows the partial barri- 
ers surrounding the largest two regions Aq and Ai (which are 
quantum mechanically accessible.) They are all constructed 
from stable/unstable manifolds of a period 4 and a period 28 
orbit. 

Using the box-counting method ll47n one can associate a 
fractal dimension Sj to the intersection F fl Aj of the chaotic 
saddle F with each of the regions Aj. The number N^^{e) of 
occupied boxes of side length e scales like N^^{e) ~ e~^^ , 
see Fig. [T] with ^o = 1-68 and Si = 1.86. To emphasize 
the difference between such dimensions close to two, the or- 
dinate is rescaled by e^. The abscissa is chosen to be e~'^ 
which will be convenient for the comparison to the fractal 
Weyl law. The increase of dj towards two when going deeper 
into the hierarchy can be qualitatively understood by adapting 
the Kantz-Grassberger relation [48] from fully chaotic sys- 
tems. The asymptotic behavior for e approaching zero will be 
discussed later in connection to the semiclassical limit. 

Quantum classification of states: We now present the 
essential quantum effect that resonance states localize pre- 
dominantly on one of the regions Aj. The closed quan- 
tum system is described by the time-evolution operator U = 
exp{-^V{q)} exp{-^T(p)} eM-±V{q)}. The corre- 
sponding open quantum system is given by /7open = PUP, 
where P is a projector on all positions not in the absorb- 
ing regions. The resonance states ip are given by /7openV^ = 
ex.-p[—i{(p — ij/2)]ilj. Regular resonance states are predomi- 
nantly located in the regular region. Chaotic resonance states 
are predominantly located in either of the hierarchical regions 
Aj, see Fig. [T] Hence, we will call them hierarchical reso- 
nance states (of region Aj.) Such a localization of chaotic 
eigenstates on different sides of a partial barrier is well known 



for closed quantum systems f^^HssllSSillll]. Chaotic eigen- 
states localized in the hierarchical re gion of a mixed phase 
space were termed hierarchical states [DTI] . They require that 
the classical flux ^ across a partial barrier is small compared 
to Planck's constant h, i. e. $ <C h, while in the opposite 
case eigenstates would be equidistributed ignoring the partial 
barrier. Quite surprisingly, we find that this condition from 
closed systems is irrelevant for hierarchical resonance states 
in open quantum systems. In the standard map at hz = 2.9 
we have ^ ~ 1/80, and for much smaller h = 1/1000 typ- 
ical resonance states still predominantly localize in one of 
the regions Aj, as shown in Fig. [T] Even for h as small as 
h = 1/12800, see Fig. [2] (insets), this is the case. This crucial 
phenomenon for our study highlights the strong impact of the 
opening. For the transition between localization and equidis- 
tribution of resonance states in the hierarchical region of open 
systems we conjecture that it depends on the ratio of escape 
rates of neighboring regions. If the escape rates and, thus, 
the corresponding quantum decay rates of resonance states lo- 
calized on different sides of the partial barrier are sufficiently 
different, such that the distance of their resonance energies ex- 
ceeds the coupling strength, the hierarchical resonance states 
do not hybridize. In contrast, for the closed system the dis- 
tance of the corresponding real eigenenergies would be much 
smaller and lead to hybridization for h < ^. This impact of 
the opening will be studied quantitatively in the future. 

For the present study it is sufficient to observe that the great 
majority of chaotic resonance states is predominantly located 
in one of the regions Aj, allowing their classification. Nu- 
merically, we use their relative local Husimi weight in Aj (in 
the case of Aq excluding the area of the opening) and discard 
states with more than 50% Husimi weight in the regular region 
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FIG. 2. (Color online) Distributions P(7) of decay rates 7 for hier- 
archical resonance states of the standard map at Ai: = 2.9 located in 
regions Aq (right) and Ai (left) for 1/h = 12800 and corresponding 
Husimi representations of typical states (insets). Short-lived states 
(7 > 7c) are not counted in the fractal Weyl law. 



and the deep hierarchical region (Aj, j > 2). The vaUdity of 
this classification is supported by the distribution of the de- 
cay rates 7 of the corresponding resonance states, see Fig. [2] 
States which are located deeper in the hierarchy have smaller 
decay rates. In fact, the two distributions for regions Aq and 
Ai have a small overlap, only. Note, that an alternative classi- 
fication of resonance states purely based on their decay rates 
7, would fail deeper in the hierarchy, as the tree-like structure 
allows for different regions Aj having strongly overlapping 
decay rate distributions. 

Hierarchical Fractal Weyl laws: For each region Aj of the 
hierarchical phase space we now relate the number ATj of hi- 
erarchical resonance states of that region to the fractal dimen- 
sion Sj of the chaotic saddle in that region. For this we use the 
fractal Weyl law of fully chaotic systems ||lllll2|], Eq. ([T]), in- 
dividually for each region Aj . This gives our main result that 
in open systems with a mixed phase space one obtains a hi- 
erarchy of fractal Weyl laws, one for each phase- space region 
Aj, Eq. Q. We stress that this result is based on the surpris- 
ing existence of hierarchical resonance states. Note, that as a 
consequence of Eq. Q the total number of long-lived hierar- 
chical resonance states is a superposition of power laws with 
different exponents and not a single power law. We will dis- 
cuss the range of validity of the hierarchical Weyl laws and 
the semiclassical limit further below. 

To give an intuitive understanding of the hierarchical fractal 
Weyl laws, let us recall the interpretation of the fractal Weyl 
law lull , and apply it in the presence of a hierarchical phase 
space. The number of quantum states localizing on a particu- 
lar phase- space region is given by the number of Planck cells 
necessary to cover the chaotic saddle in that region. Using 
the scaling, N^^{e) ~ e~^^ , of the number of boxes N^^ to 
cover the chaotic saddle in region Aj and the identification 
of the box area e'^ with the Planck cell area h directly leads 
to Eq. (|2]). In analogy to fully chaotic open systems ll49ll . 
the fine- scale structure of the hierarchical resonance states is 
dominated by the backward trapped set for right eigenstates 
of ^open and the forward trapped set for left eigenstates, see 
insets in Fig. O 

The numerical investigation of the standard map supports 
the existence of hierarchical fractal Weyl laws, as we now 
show. By the classification of resonance states we are able 
to determine the number J\fj {h) of long-lived hierarchical res- 
onance states associated to a particular region Aj depending 
on h. We restrict ourselves to the consideration of small h 
such that ^/h > 10 where quantum mechanics can very 
well mimic classical transport in phase space [41]. Short- 
lived states are discarded by defining an arbitrary cut-off rate 
7c = 1, as usual for the fractal Weyl law [12]. Here this af- 
fects resonance states of the outermost region Aq only. We 
obtain distinct behavior for each counting function JVj, see 
Fig. [T] Here the ordinate is rescaled by h and the abscissa 
is chosen to be h~^, corresponding to the previous classical 
rescalings. For small h one observes the power-law scaling 
of Eq. ([2| and good agreement with the box-counting results 
for the fractal dimensions Sj of F fl Aj. Note that we fit- 
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FIG. 3. (Color online) Hierarchical fractal Weyl laws J\fj vs. h~^ 
(filled symbols) counting hierarchical resonance states in the outer 
(Ao, triangles), central (Ai, circles), and inner (A2, squares) chaotic 
regions for the hierarchical model system, and corresponding typical 
Husimi representations (right) f or /i = 1/1115. They are compared 
to the box-counting scaling N^^ vs. e~'^ (open symbols) giving dis- 
tinct fractal dimensions 6j of the corresponding regions of the chaotic 
saddle (inset). Here we used /o = 1.75, /i = 1.55, and /2 = 0.8 
which again are of order one. 



ted prefactors fj to the quantum results to better demonstrate 
that they scale with power-law exponents in agreement with 
the classical counterparts. Both prefactors fj are of order one 
(/o = 2.6, /i = 0.85). Figure [T] confirms for two regions Aj 
of the standard map that they give rise to hierarchical fractal 
Weyl laws. 

Hierarchical model system: To verify the hierarchical frac- 
tal Weyl laws for more than two regions, we suggest the fol- 
lowing system that models the hierarchical structure of partial 
barriers in a generic mixed phase space, similar in spirit to 
a one-dimensional model |29] and a Markov chain [26]. In 
contrast to the standard map, this designed model quantum 
mechanically allows for studying three regions numerically. 

We first define a composed symplectic map C o M on the 
phase space F = [0, 1) x [0, 1). It models b partial barri- 
ers at the positions qi < • • • < qb sls straight lines in p- 
direction, giving a decomposition of F into b -\- 1 regions 
Aj = [qj.qj^i) X [0, 1) with qo = and q^^i = 1. The map 
M describes the uncoupled dynamics being sufficiently mix- 
ing in each region Aj . We choose the standard map at kicking 
strength k, = 10 acting individually on each of the regions Aj 
after appropriate rescaling. The map C couples these regions 
mimicking the turnstile mechanism of a partial barrier with 
flux ^j by exchanging the areas [qj — ^j , g^ ) x [0, 1) with their 
neighboring areas [qj^qj + ^j) x [0, 1). Finally, we open the 
system by defining the absorbing region [0, ^0) x [0, 1) and 
the corresponding escape map O, leading to the hierarchical 



model map O oC o M oO. The barrier positions qj are cho- 
sen such that the area fJ^{Aj) = qj^i — Qj of the regions Aj 
obeys the scaHng jj.{Aj)/fi{Aj_^i) = a with some a > 1. We 
choose the scaHng behavior of the fluxes as ^j/^j-^i = (j) 
with some (j) > a, such that even for an arbitrary deep hierar- 
chy, the size of the coupHng regions wiU not exceed the size 
of the regions Aj. In this paper we use the parameters 6 = 2, 
Of = 2, ^ = 4, and ^o//i{Ao) = 1/4. Hence, we have the 
fluxes $1 = 1/28 and <l>2 = 1/112. 

Figure [3] shows the results for the hierarchical model sys- 
tem: We obtain the fractal dimensions So = 1.69, Si = 1.94, 
and S2 = 1.99. Quantum mechanically, we again find hier- 
archical resonance states predominantly localizing on one of 
the regions Aj, even though h <C ^1, ^2- There number fol- 
lows the proposed hierarchical fractal Weyl laws according to 
Eq. ©. 

Semiclassical limit: We now discuss the range of validity 
of the individual fractal Weyl law, Eq. (|2]), for some region 
Aj as a function of decreasing Planck's constant h. Once h is 
small enough to resolve this region Aj, i.e. h ^ Aj, there are 
four relevant regimes: (i) For h larger than the greatest flux ^j 
across its surrounding partial barriers, h > ^j, one has reso- 
nance states localized on region Aj with just a small coupling 
to other regions, as for closed systems ||33l|4l|]. Consequently, 
the number of resonance states in this regime scales with the 
usual Weyl law as h~^. The resonance states do not resolve 
the fine-scale structure of the chaotic saddle, (ii) For h smaller 
than the flux, h <^j, the resonance states still localize in re- 
gion Aj . Furthermore, they begin to resolve the fractal struc- 
ture of the chaotic saddle. This is the main regime discussed in 
this paper and described by Eq. ^ with a fractal dimension 
Sj of the intersection of the chaotic saddle with region Aj . 
(iii) For much smaller values of h, we expect that hierarchi- 
cal resonance states start to extend over several regions Aj . 
This expectation is based on the fact that the effective frac- 
tal dimensions of these regions coincide for sufficiently small 
side length e of the boxes, due to their coupling across par- 
tial barriers with the largest dimension dominating [29 J. As a 
consequence we expect that the quantum decay rates become 
identical leading to a hybridization of the resonance states. 
This classical scale e implies the beginning of this regime on 
the quantum scale h = e^. Owing to an extremely slow e 
dependence of the effective fractal dimensions these quantum 
mechanical consequences are beyond current numerical inves- 
tigations, (iv) In the semiclassical limit, h -^ 0, we presume 
an overall Weyl law for the hierarchical region with the num- 
ber of resonance states scaling as h~^. This is based on the 
fact that asymptotically (e -^ 0), the fractal dimension is the 
same in all regions Aj [29], and that for an infinite hierarchi- 
cal structure of partial barriers, the fractal dimension of the 
chaotic saddle equals two lisoll . 

Discussion: In this paper we have shown that the hierarchi- 
cal fractal Weyl laws, Eq. Q, describe the important regime 
(ii) where hierarchical resonance states predominantly local- 
ize on one of the regions Aj and resolve the fractal structure of 
the chaotic saddle. Note that Eq. Q also applies to the other 



regimes by choosing the phase-space regions according to the 
predominant localization of resonance states. 

There is also a related Weyl law with fractional exponent 
in the presence of a hierarchical phase space even for closed 
systems. In Ref. [I37l1 it was shown that the number of hierar- 
chical states scales Sish^~^^^, relating this Weyl law to the ex- 
ponent 7 of the power-law decay of the Poincare recurrence- 
time distribution P{t) ~ t~^ with 7 > 1. These hierarchical 
states localize on regions Aj with flux ^j < h. In the open 
system they correspond to resonance states of all regions in 
regime (i), which, thus, fulfill the same fractional Weyl law 
~ /i,i-i/7 xhe idea of connecting the fractal Weyl law with 
the power-law decay of P{t) was also studied in Ref. ll44ll . 
There a fractal Weyl law ^ h^~^ was found for the number of 
resonance states fulfilling the cut-off criterion 7 < 7c ~ /i. It 
can be shown that this includes resonance states from regimes 
(i) and (ii) when applied to a hierarchical phase-space struc- 
ture. In contrast, the hierarchical fractal Weyl laws, Eq. (O, 
focus on resonance states of regime (ii). There, even though 
h < ^j, the resonance states predominantly localize on one 
region Aj. Furthermore, these resonance states resolve the 
fine-scale structure of the chaotic saddle. A future challenge 
is the study of fractal Weyl laws in higher dimensional sys- 
tems with a generic phase space. 
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